Measurement of Parton Distributions of Strange Quarks in the Nucleon from 
Charged-Kaon Production in Deep-Inelastic Scattering on the Deuteron 
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The momentum and helicity density distributions of the strange quark sea in the nucleon are 
obtained in leading order from charged-kaon production in deep-inelastic scattering on the deuteron. 
The distributions are extracted from spin-averaged K multiplicities, and from K ± and inclusive 
double-spin asymmetries for scattering of polarized positrons by a polarized deuterium target. The 
shape of the momentum distribution is softer than that of the average of the u and d quarks. In the 
region of measurement 0.02 < x < 0.6 and Q 2 >1.0 GeV 2 , the helicity distribution is zero within 
experimental uncertainties. 

PACS numbers: 13.60.-r, 13.88.+e, 14.20.Dh, 14.65.-q 



Parton distribution functions (PDFs) form the basis 
for the description of the flavor structure of the nucleon. 
The spin-averaged parton distribution functions q(x) of 
quarks and antiquarks of flavors q = (u,d,s) [3,11] de- 
scribe the quark momentum contributions, where x is 
the dimcnsionlcss Bjorkcn scaling variable representing 
the momentum fraction of the target carried by the par- 
ton in a frame where the target has "infinite" longitudinal 
momentum. They are sums of the number densities of 
the quarks q^(x) [q^±(x)} with the same [opposite] he- 
licity as that of the nucleon. The differences, or helicity 
distributions, Aq(x) = q^(x) — q^(x) describe the flavor 
dependent contributions of the quark spins to the spin 
of the nucleon. The features of the parton distributions 
reflect the QCD dynamics of the constituents. Because 
strange quarks are objects which reflect directly proper- 
ties of the nucleon sea, they arc of special interest. Their 
distributions are also important because of their impact 
on quantitative calculations of certain key short-distance 
processes at hadron colliders, and their implications for 
the measurement of the Weinberg angle in deep- inelastic 
scattering (DIS) of neutrinos @, S]- 

In the absence of significant experimental constraints, 
current global QCD fits of PDFs [5, 6] assume the strange 
quark and antiquark momentum distributions s(x) and 
s(x) to be given by s(x) = s(x) = r\u(x) + d(x)]/2 
with r w 1/2 at some low factorization scale. Mea- 
surements of neutrino and antineutrino production of 
dimuons @, 1, 0, M, EIJI1 EE Q 01 provide useful 
but limited information [161 ] on the normalization and 
shape of the distribution s(x) +s(x). In these experi- 
ments, extraction of the strange quark distributions re- 
quires knowledge of the charm quark mass, the charm 
hadron semileptonic branching ratio, and the "Peterson 
fragmentation parameter" [l7| that describes the kine- 
matic dependence of the charm fragmentation function. 
These quantities together with the strange parton distri- 
butions themselves are fitted simultaneously in the ex- 
traction procedure. Much of the information on proper- 
ties of the helicity distribution of strange quarks is based 
on the analysis of inclusive DIS and hypcron decay under 
the assumption of SU(3) symmetry among the structures 



of the octet baryons. In these inclusive experiments (l8j 
the first moment of the helicity distribution for strange 
quarks is one of the principal results. The most precise 
recent value is -0.103 ± 0.007(exp.) ± 0M3(theor.) ± 
0.008(euoZ.) in LO [ljj. A full 5-flavor decomposition 
using Hermes semi-inclusive DIS [2^] data from proton 
and deuteron targets, although not sensitive to As(x), 
yielded As = 0.028 ± 0.033 ± 0.009 for the first partial 
moment of the strange quark helicity density in the mea- 
sured range 0.023 < x < 0.3. A separate "isoscalar" ex- 
traction of As + As from DIS data on the deuteron alone 
gave As + As = 0.129 ± 0.042 ± 0.129 in the measured 
range where the large systematic uncertainty reflected 
lack of knowledge of kaon fragmentation functions. 

This letter reports a new isoscalar extraction of s(x) + 
s(x) and A(s(x) +s(x)) based on the same Hermes data 
obtained from polarized DIS on a deuterium target. The 
measurement reported here is complementary to the neu- 
trino results, and is the first extraction of s(x) + s(x) 
in charged lepton DIS. Because strange quarks carry no 
isospin, the strange seas in the proton and neutron can be 
assumed to be identical. In the deuteron, an isoscalar tar- 
get, the fragmentation process in DIS can be described by 
fragmentation functions that have no isospin dependence. 
Aside from isospin symmetry between proton and neu- 
tron, the only symmetry assumed is charge-conjugation 
invariance in fragmentation. For the isoscalar deuteron 
in Leading Order (LO), the inclusive unpolarized (U) 
electron scattering cross section in terms of the parton 
distributions Q(x) = u(x) + u(x) + d(x) + d(x) and 
S(x) = s(x) +s(x) takes the form 



d 2 N DIS {x) 



ICu(x,Q 2 )[5Q{x)+2S(x)} , (1) 



dxdQ 2 

where )Cu{x,Q 2 ) is a kinematic factor containing the 
hard scattering cross section. The weak logarithmic 
dependence of the PDFs on — Q 2 , the squared four- 
momentum of the exchanged virtual photon, has been 
suppressed for simplicity. Applying the same LO formal- 
ism to the semi-inclusive cross section for charged kaon 
production, irrespective of charge, hereafter designated 
as K gives 

d 2 N K {x) 



dxdQ 2 
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Q(x) I V% (z)dz + S(x) I £>£ (z)dz 



vK", 



(2) 
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where z = Ehjv with v and Eh the energies of the vir- 
tual photon and of the detected hadron in the target rest 
frame, V* (z) = 4D^(z)+D^(z) and V§ (z) = 2D* (*). 
The fragmentation function D^{z) describing the num- 
ber density of charged kaons from a struck quark of flavor 
q is integrated over the measured range of z. Combining 
Eqs. (|H2p and neglecting the term 2S(x) compared to 
5Q(x), it follows immediately that 



S(x) I Pf \z)dz ~ Q(x) 



d 2 N K (x) 
d 2 N DIS {x) 



V$(z)dz 



(3) 

Eq. [3] is the basis for the extraction of the quantity 
S(x)fD§(z)dz. 

The data were recorded with a longitudinally nuclear- 
polarized deuteron gas target internal to the E = 27.6 
GeV Hera positron storage ring at Desy. The self- 
induced beam polarization was measured continuously 
with Compton backscattering of circularly polarized laser 
beams [23, [22| . The open-ended target cell was fed by an 
atomic-beam source based on Stern-Gerlach separation 
with hyperfine transitions. The nuclear polarization of 
the atoms was flipped at 90s time intervals, while both 
this polarization and the atomic fraction inside the target 
cell were continuously measured The average value 
of the deuteron polarization was 0.845 with a fractional 
systematic uncertainty of 3.5%. 

Scattered beam leptons and coincident hadrons were 
detected by the Hermes spectrometer j24j. Leptons 
were identified with an efficiency exceeding 98% and a 
hadron contamination of less than 1% using an elec- 
tromagnetic calorimeter, a transition-radiation detec- 
tor, a preshower scintillation counter and a ring-imaging 
Cerenkov (RICH) detector [25j]. The dual-radiator 
RICH was also used to identify charged kaons. Events 
were selected subject to the kinematic requirements 
Q 2 > lGcV 2 , W 2 > lOGeV 2 and y < 0.85, where 
W is the invariant mass of the photon-nucleon system, 
and y = v/E. Coincident hadrons were accepted if 
0.2 < z < 0.8 and x F « 2p L /W > 0.1, where p L is the 
longitudinal momentum of the hadron with respect to the 
virtual photon direction in the photon-nucleon center of 
mass frame. The Bjorken x range of measurement was 
0.02-0.6. 

The charged kaon multiplicity was extracted by sum- 
ming over the kaon yields for the two beam-target po- 
larization states. An event weighting procedure was 
used to correct for RICH kaon identification inefficien- 
cies. The effects of QED radiation, instrumental reso- 
lution, and acceptance were simulated [II Hill], and 
corrections were applied to the data for each polarization 
state usin g a technique that unfolds kinematic migration 
of events [19( . The results are presented in Fig. [1] The 
trends in the data were not reproduced (see dotted curve 
in Fig. [1]) by fitting the points using the Cteq6l S 
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FIG. 1: The multiplicity corrected to 4tt of charged kaons in 
semi-inclusive DIS from a deuterium target, as a function of 
Bjorken x. The continuous curve is calculated from the curve 
in Fig. fusing Eq. [3j The dashed(dash-dotted) curve is the 
nonstrange(strange) quark contribution to the multiplicity for 
this fit. The dotted curve is the best fit to f T>§ (z)dz using 
CteqGl PDFs. The error bars are statistical. The band 
represents the systematic uncertainties. The values of (Q 2 ) 
for each x bin are shown in the lower panel. 



and J T>g (z)dz as free parameters. In view of the paucity 
of reliable data on S(x), it was assumed instead that it 
is unknown, and the analysis was carried out extract- 
ing the product S(x) J Dg(z)dz in LO. For x > 0.15 
the multiplicity is constant at a value of about 0.080, 
implying that S(x)/Q(x) is constant. For this analy- 
sis S(x) is assumed to be negligible at large x from 
which it follows that S(x) = for x > 0.15 and that 
Jo 2 ^q( z )^ z = 0.398 ± 0.010, in excellent agreement 
with the value 0.435 ±0.044 obtained for Q 2 = 2.5 GeV 2 
from the most recent global analysis of fragmentation 
functions [3(| ■ The value 0.398 was then used in Eq. ([3]) 
together with values of Q{x) from Cteq6l and the mea- 
sured multiplicities to obtain the product S{x) J Vg{z)dz 
shown in Fig. O A small iterative correction was made 
to account for the neglect of the 2S(x) term in Eq. [T] 
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FIG. 2: The strange fragmentation product 

S(x,Q 2 )fT>s(z)dz obtained from the measured Hermes 
multiplicity for charged kaons at the (Q 2 ) for each bin. The 



curve is a least squares fit of the form x 



-0.863 -x/0. 0487 



(1-x). 



strange quark PDFs in Eqs. ffl and El with / V%(z)dz The band re P resen ts systematic uncertainties. 
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FIG. 3: The strange parton distribution xS(x) from the 
measured Hermes multiplicity for charged kaons evolved to 
Qo = 2.5 GeV 2 assuming / £># (z)dz = 1.27±0.13. The solid 



curve is a 3-parameter fit for S(x) 



-0.924 -x/0. 0404 



(1-x), 



the dashed curve gives xS(x) from Cteq6l, and the dot-dash 
curve is the sum of light antiquarks from Cteq6l. 



The result for the product together with a fit of the form 
x~ ai e~ x / a2 (1 — x) is shown in Fig. O and leads to the 
continuous curve in Fig. [T] 

The improved fit (continuous curve in Fig. [T]) to the 
multiplicity is an indication that the actual distribution 
of S(x) is substantially different from the average of those 
of the nonstrange antiquarks. To explore this point, the 
Hermes result for S(x)J T>g(z)dz has been evolved to 
Qq = 2.5 GeV 2 . The Q 2 evolution factors were taken 
from Cteq6l and the fragmentation function compila- 
tion given in [30| . Consideration of corrections to the 
evolution due to higher twist contributions is not nec- 
essary, since higher twist effects are expected to be sig- 
nificant [3~J ] only for larger values of x where the ex- 
tracted distribution of xS(x) vanishes. The distribution 
of xS(x) was obtained from S(x) J T>g (z)dz by dividing 
by / V§(z)dz = 1.27±0.13, the value at Q 2 = 2.5 GeV 2 
given in [30j] . The results are presented in Fig. [3l The 
normalization of the Hermes points is determined by the 
value of J T>g (z)dz assumed. However, whatever the nor- 
malization, the shape of xS(x) implied by the Hermes 
data is incompatible with xS (x) from Cteq6l as well as 
the assumption of an average of an isoscalar nonstrange 
sea. The absence of strength above i w 0.1 is clearly dis- 
crepant with Cteq6l, while deviations from the Cteq6l 
prediction at low x could be, in part, a manifestation of 
higher order processes. 

In the isoscalar extraction of the helicity distribution 
AS(x) = As(i)+As(i), only the double-spin asymmetry 
A^ d (x, Q 2 ) for all charged kaons, irrespective of charge, 
and the inclusive asymmetry A\\ d (x,Q 2 ) are used. In 
LO, the inclusive and the charged kaon double-spin(LL) 
asymmetries are determined by the relations 



d 2 N DIS (x) 



MW dxdQ 2 

JC LL (x,Q 2 )[5AQ(x) 
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FIG. 4: Lepton-nucleon polarized cross section asymmetries 
Ay d for inclusive DIS and A^ d for semi-inclusive DIS by 
a deuteron target as a function of Bjorken x, for identified 
charged kaons. The error bars are statistical, and the bands 
at the bottom represent the systematic uncertainties. 



where K-ll is a kinematic factor, and 
d 2 N K (x) 



ICll(x,Q z )x 



dxdQ 2 

AQ(x) J V^{z)dz + AS(x) J V 1 §{z)dz 



(5) 



2AS(x)\ 



(4) 



Eqs. (|4I5[) permit the simultaneous extraction of the he- 
licity distribution AQ(x) = Au(x) + Au(x) + Ad(x) + 
Ad(x) and the strange helicity distribution AS(x) = 
As(x) + As(i). The nonstrange integrated fragmenta- 
tion function needed for a LO extraction of AS(x) was 
extracted from the multiplicity analysis of the same data. 

The semi-inclusive asymmetries A^ d were derived from 
the kaon spectra measured for each target polarization. 
The target polarization was corrected for the D-wave ad- 
mixture in the deuteron wave function by applying the 
correction term (1 — 1.5u;rj) where u> = 0.05 ± 0.01 [32I ]. 
The corrected asymmetries are shown in Fig. 01 The in- 
clusive asymmetries A^^x) were corrected for effects of 
QED radiation and instrumental smearing with the same 
procedures described above for the spin dependent kaon 
multiplicities. Contributions to the systematic uncertain- 
ties in the asymmetries include those from the beam and 
target polarizations, and the ne glec t of the transverse 
spin structure function 52(2;) ~ [331 ] . and for A^ d from 
those of RICH kaon identification. 

The quark helicity distributions were extracted from 
the measured spin asymmetries and A^ d (x) in an 

analysis based on Eqs. (|4l5p . The value of J T>g (z)dz = 
1.27 ± 0.13 was used to extract AS(x). The results are 
presented in Fig. [5] The strange helicity distribution also 
agrees well with the less precise results of and is 
consistent with zero over the measured range. 

The first moments of the helicity densities in the 
measured region are presented in Tab. |TJ The result 
for AQ over the measured range is consistent with the 
value 0.381 ± 0.010(stat.) ± 0.027(sys.) for the full mo- 
ment previously extracted from Hermes gi t d data 
The value of AS measured here is not in serious dis- 
agreement with -0.0435 ± 0.010(stat.) ± 0.004(sys.) ex- 
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FIG. 5: Nonstrange and strange quark helicity distributions 
at Ql = 2.5 GeV 2 , as a function of Bjorken x. The error 
bars are statistical, and the bands at the bottom represent 
the systematic uncertainties. The curves are the LO results 
of Leader et al. [3^ from their analysis of world data. 



tracted from the inclusive Hermes measurements. The 
value for the partial moment of the octet combination 
Aq$(x) = AQ(x) - 2AS(x), included in Tab. [fl is sub- 
stantially less than the value of the axial charge as = 
Ao 8 = Jg 1 Aq 8 (x)dx = 0.586 ± 0.031 extracted from 
the hyperon decay constants by assuming SU(3) sym- 
metry [34J . Possible explanations for the deficit observed 
for Aqg, include violation of SU(3) symmetry or missing 
octet strength at values of x below the measured range. 
The substantial deviation observed in the shape of S(x) 
from that of the light sea quarks is a clear manifestation 
of violation of SU(3) symmetry 35|, 3f| 37 1 in the strange 
quark sector. 

In conclusion, inclusive and semi-inclusivc-charged- 
kaon spin asymmetries for a longitudinally polarized 
deuteron target have been analyzed to extract the LO 
parton distributions of the strange sea in the proton. The 
partial moment of the nonstrange fragmentation func- 
tion needed for the LO analysis has been extracted di- 
rectly from the same data. The values for the PDFs pre- 
sented in this paper arc available at the Hermes web site 
( |http: / / www-hermes.desy.de| . The momentum densities 
are softer than previously assumed. The helicity densi- 
ties are consistent with zero and the partial moment of 
the octet axial combination is observed to be substan- 
tially less than the axial charge extracted from hyperon 
decays under the assumption of SU(3) symmetry. 



TABLE I: First moments of various helicity distributions in 
the Bjorken x range 0.02-0.6 at a scale of Q 2 , = 2.5 GeV 2 . 
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Moments in measured range 


AQ 


0.359 ± 0.026(stat.) ± 0.018(sys.) 


AS 


0.037 ± 0.019(stat.) ± 0.027(sys.) 


Aqg 


0.285 ± 0.046(stat.) ± 0.057(sys.) 
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